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ABSTRACT 


The  statistical  properties  of  rain  attenuation  which  dominate 
propagation  at  frequencies  above  10  GHz  are  required  for  planning 
radao  systems  at  these  frequencies.  These  propoerties  are  studies  for 
three  frequencies,  11.2,  12.7  and  18. 76  GHz  and  various  link  lengths, 
2.5,  5»  7»5j  10,  12.5  and  15  mi  at  six  locations  in  Canada.  Two 
approaches  were  used  to  simulate  the  attenuation;  one  utilizes 
weather  radar  data  in  the  AZLOR  format,  a rectangular  coordinate 
variation  of  PPI;  the  other  utilizes  rainguage  records  and  radiosonde 
data  and  the  concept  of  a sj*nthetic  stom.  Rain  attenuation  statistics 
are  developed  from  these  two  approaches  to  describe  the  probalbility  of 
occurrence  of  an  attenuation  event  for  single  links  of  various  length, 
pairs  of  parallel  links  with  various  separation  distance  and  for  pairs 
of  links  on  the  same  routine. 


PREFACE 


The  aim  of  this  report  is  to  collect  in  one  place  and  in  printed 
form  the  results  of  three  years  of  work  on  rain  attenuation  statistics. 
Nothing,  or  close  to  nothing,  has  been  done  here  toward  the  explanation 
or  rational  interpretation,  of  these  results  in  terms  of  storm  modelling 
or  climatology,  although  these  results  contain  much  information  on  these 
topics.  These  results  shoiild  just  be  consiedered  for  what  they  are: 
observations  of  the  precipitation  process  translated  in  a form  suitable 
for  the  planning  of  radio  systems  across  Canada. 

I am  grateful  to  the  Atmospheric  Environment  Service  for  rain- 
guage  and  upper-wind  data,  and  to  McGill  Radar  Weather  Observatory  for 
radar  records,  (The  basic  radar  there  is  on  loan  from  the  Air  Force 
Cambridge  Research  Laboratories  under  their  contract  F19628-73-C-0052.) 

I wish  to  acknowledge  the  contribution  of  Hr,  Chul-Un  Ro  v/ho 
carried  on  the  analysis  of  the  rainguage  data  for  Montreal  and  the 
other  locations  in  Canada  during  his  work  toward  a Master’s  degree  in 
Meteorology. 

I am  deeply  indebted  to  Prof,  J.S,  Marshall  and  to  the  other 
members  of  the  Stormy  VJeather  Group;  in  particular  to  ray  friends 
Dr.  R.R,  Rogers  and  I,  Za\\radzki  wlio  had  many  suggestions  and  useful 
criticisms.  Dr.  K.S,  McCormick  and  D.S,  Strickland  of  CRC  (which 
supported  the  work)  have  been  most  helpful  and  understanding. 

Mr. Daniel  Salomon  did  all  the  computer  work  and  his  contri- 
bution is  greatly  appreciated. 
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RAIN  attenuation  STUDIES 


1.  INTRODUCTION 

For  the  planning  of  radio  systems 
making  use  of  frequencies  above  10  CHz  the 
knowledge  of  the  statistical  properties  of 
rain  attenuation  which  dominate  the  propa- 
gation at  these  frequencies  Is  required. 

The  basic  parameter  In  the  design  of 
a radio  system  Is  the  outage  time  of  the  system 
which  can  be  due  to  various  causes,  among  which 
Is  attenuation  by  rain.  The  outage  time  due 
to  rain  attenuation  depends  on  the  possibility 
of  the  system  to  recover  from  the  noise  on 
attenuated  signal.  This  possibility  Is  limited 
to  a given  threshold,  which  Is  a function  of 
the  system  characteristics.  The  outage  time 
Is  then  the  total  time  for  which  the  rain 
attenuation  exceeds  this  threshold. 

This  work  attempts  to  describe  the 
statistical  properties  of  the  time  for  which 
rain  attenuation  exceeds  a given  value  as  a 
function  of  the  frequency  and  the  hop  length. 
These  properties  are  studied  for  three  fre- 
quencies (11. 2j  12. 7j  18.76  GHz)  and  various 
hop  lengths  (2.5j  5j  7.5j  10;  12. 5>  15  I^i)* 

Results  are  given  for  six  locations  across 
Canada.  They  are:  Edmonton,  Winnipeg, 

Thunder  Bay,  Toronto,  Montreal  and  Halifax. 

Two  approaches  have  been  used,  the 
first,  described  In  chapter  2,  makes  use  of 
the  McGill  weather  radar  data  for  1971,  jointly 
with  the  data  of  an  experimental  link  at  11.2 
GHz  installed  by  BNR  In  Ormstown  (fig.  l).  In 
this  experiment  the  radar  data  of  the  summer 
19~1  have  been  used  to  simulate  the  attenuation 
process  over  a dense  network  of  links  of  various 
lengths,  and  the  observed  statistics  from  the 
BNTt  link  have  been  used  to  calibrate  the  simu- 
lation. This  experiment  has  given  unique  Infor- 
mation about  extreme  events,  and  thanks  to  the 
very  large  sample  size  of  the  data,  events 
with  probability  of  10"^  have  been  evaluated. 
Further  it  has  provided  unique  information 
about  joint  probabilities  of  attenuation  for 
parallel  links  such  as  required  tor  the  design 
of  twin  route  systems. 

A second  approach,  described  in  chap- 
t-  r consists  in  the  develop-ipnt  of  a methol 
"i  ->L"ulatine  the  attenuation  prf.ccss  from 


raingauge  records  and  radiosonde  data.  The 
method  Is  based  on  the  concept  of  a synthetic 
storm  first  Introduced  by  Hamilton  and  Marshall. 
The  method  has  been  tested  with  the  aid  of  the 
McGill  Weather  Radar  and  the  BNR  link  and  then 
applied  to  ten  years  of  data  for  the  six  loca- 
tions across  Canada  (chapter  4).  The  results 
allow  study  of  the  year  to  year  variations  of 
the  statistics  of  interest  and  will  be  of  help 
in  planning  radio  systems  across  Canada. 


2.  ATTENUATION  BY  RADAR 


2.1  Introduction 

The  purpose  of  this  work  has  been  to 
obtain  attenuation  statistics  for  the  Montreal 
region.  In  order  to  design  a radio  system  It 
is  desirable  to  know  the  probability  of  occur- 
rence of  attenuation  events  for  one  link  down 
to  a very  low  fraction  of  time  (typically  10"® 
of  one  year).  For  the  Montreal  region,  one 
microwave  link  operating  for  a one  year  period 
provides  a sample  of  sufficient  size  to  con- 
tain cvenp  with  a typical  probability  of 
about  1Q"°.  In  a first  approximation  it  would 
be  necessary  to  operate  some  100  links  a year 
to  reach  the  desired  probability  level,  and 
this  is  just  for  one  frequency  and  one  hop 
length.  It  is  obvious  that  this  approach  is 
not  practical  or  economical.  As  an  alternative, 
the  information  requested  can  be  obtained  empir- 
ically from  several  years  of  raingauge  records 
(chapter  3)  or  from  radar  data.  A radar  can 
cover  a large  geographical  region,  thus  provid- 
ing a large  number  of  paths  over  which  atten- 
uation events  can  be  calculated.  In  this  work 
this  approach  has  been  carried  out,  and  atten- 
uation statistics  have  been  evaluated  at  11.2, 
12.7,  15.7  GHz.  These  statistics  are  for 
single  links  of  various  length,  for  pairs  of 
parallel  links  with  various  separation  dis- 
tances and  for  a pair  of  links  on  the  sam»> 
route.  With  the  method  oinployecl  no  inforria- 
tion  about  the  duration  of  the  events  is  ob- 
tained, although  this  inf  or;  .a  C ion , from  the 
raingauge  data,  is  given  in  ch->'ter  Radar 


simulation  of  attenuation  has  been  used  and 
studied  by  various  authors  (l»2,3i)»  There 
are  two  major  problems  that  arise  In  this  kind 
of  approach:  calibration  of  the  radar  and 

evaluation  of  absolute  probabilities.  They 
will  be  discussed  In  2.3  snd  2.4« 

The  data  used  In  this  study  have  been 
provided  by  the  McGill  Weather  Radar  and  they 
cover  the  summer  months  of  1971  • A brief 
description  of  the  radar  follows. 


2.2  The  McGill  Weather  Radar 

During  1971  the  main  features  of  the 
radar  were: 

Frequency  2,880  MHz 

Antenna  Diameter  30*foot  paraboloid 

(0,  / = 0.8°) 

Peak  Power  2 MW 

Pulse  Repetition  Frequency  60  Hz 

Pulse  length  0.3  sec  (pulse  com- 
pressing a 4 >1  sec  pulse  with  a 
factor  13) 

Sensitivity  - 104  dBm  (equivalent  to 
a rainfall  rate  of  0.2  mn/hr  at 
200  km) 

Dynamic  Range  80  dB. 

Every  S minutes  the  radar  scans  a 
hemisphere  with  a radius  of  200  km.  The 
scanning  consists  of  successive  azimuth  rota- 
tions at  different  elevation  angles  from  0.3° 
to  32° • Each  rotation  Is  completed  In  10 
seconds.  The  measured  reflectivity  Is  stored 
on  35  mm  films  with  3 dB  thresholds.  Use  Is 
made  of  the  range  resolution  of  30'meters  In 
combining  the  data  from  several  contiguous 
contributing  volumes  (3  at  the  closest  range 
of  25  km  and  80  at  the  furthest  range  of 
200  km)  and  then  peak  detecting  In  order  to 
reduce  random-phase  fluctuations  (4»3i6,7)» 

The  data  format,  AZLOR,  Is  on  rectan- 
gular coordinates,  azimuth  (600  values)  on 
abscissa  and  logarithm  of  range  on  ordinate 
(200  values  going  from  23  to  200  km)  (S). 

Although  this  radar  system  has  not 
been  designed  for  attenuation  studies.  It 
Is  very  adequate  for  the  purpose  and  some 
of  its  features  are  very  helpful  In  this 
respect. 

The  AZLOR  data  format,  although  It 
distorts  shapes  and  does  not  conserve  areas. 

Is  very  rational  and  is  well  suited  for  the 
attenuation  simulation  carried  on  by  computer. 
The  main  advantage  over  PPI  is  that  the  data 
are  relative  to  square  areas  so  that  links 
can  be  simulated  in  radial  and  tangential 
direction  easily  and  promptly. 


The  antenna  program  limits  the 
sampling  time  to  5 minutes,  when  one  Is 
Interested  in  a simulation  over  a place  as 
close  to  the  ground  as  possible,  but  offers 
the  possibility  of  using  a plane  at  higher 
elevation  In  the  frequent  case  when  the 
weather  at  lower  elevation  is  masked  by 
ground  clutter,  thus  Increasing  substantially 
the  sample  size.  This  procedure  must  be 
applied  with  some  caution  in  relation  to  the 
altitude  of  the  echo,  possible  presence  of 
bright  band  and  generally  altitude  of  the 
freezing  level,  to  avoid  a contamination  of 
the  data  with  precipitation  other  than  rain 
for  which  attenuation  properties  are  not  yet 
well  known. 


2.3  Data  collection  and  analysis 

The  1971  summer  data  1>ave  been  used 
to  calculate  attenuation  statistics.  For 
this  purpose,  out  of  the  elevation  cycles  of 
the  antenna,  one  arimuth  scan  was  chosen  with 
as  low  an  elevation  angle  as  possible  In 
order  to  consider  echoes  below  freezing  level 
but  high  enough  to  avoid  ground  clutter  ef- 
fects. Except  for  a few  cases  the  elevation 
angle  used  was  0.8°  above  the  horizontal.  In 
order  to  reduce  Che  amount  of  data  to  be 
analyzed,  a threshold  was  applied  to  the 
original  reflectivity  data  such  that  only 
attenuation  events  over  10  dB  were  considered. 
For  each  scan  the  reflectivity  data  were  read 
from  the  film  and  punched  on  cards.  A pro- 
gram made  corrections  for  peak  reading  as  a 
function  of  Che  number  of  independent  samples, 
which  depends  on  range.  The  program  then 
transformed  the  data  from  reflectivity  Co 
attenuation  using  the  following  relations: 

Y(db/mi)  «=  .0003795  Z’’^^^(mmVn>^)  @ H*2  GHz 

Y = .0006431  @ 12.7  GHz 

Y = .00258  @ 18. 7 GHz 

(2.3.1) 

The  data  were  integrated  over  a dense 
grid  of  paths  of  various  lengths  and  separation 
distances.  The  attenuation  events  were  weighted 
in  function  of  range  because  of  the  logarithmic 
distribution  in  range  and  then  they  were  collec- 
ted In  histograms. 

The  network  was  composed  of  parallel 
and  perpendicular  links  directed  along  and 
across  range.  The  network  of  links  did  not 
have  a fixed  geographical  location  but  it  was 
simply  superimposed  on  a patch  of  weuher  so 
that  the  maximum  possible  number  of  attenuation 
events  was  observed. 
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In  this  way  the  possibility  of  obtain- 
ing absolute  probabilities  is  apparently  lost; 
but  this  shouldn’t  be  cause  of  grief  because 
this  possibility  has  been  lost  already  for 
more  helpless  reasons  than  inadequate  way  of 
processing  the  data*  In  fact  the  radar  wasn’t 
working  all  the  timei  for  various  reasons: 
power  failures,  equipment  failures,  maintenance 
requirements,  among  which  one  in  particular  is 
annoying,  i.e*  the  compulsion  technicians  have 
to  run  calibration  patterns  in  the  middle  of 
the  most  Interesting  convective  activity  and 
then  to  run  them  again  after  a short  time. 

Over  the  whole  summer,  the  usable  data  for 
this  analysis  covered  six  different  days  for 
a total  of  22  hours  and  20  minutes,  which  is 
a small  part  of  the  whole  summer. 

One  must  bear  in  mind  that  the  data 
collected  were  subject  to  the  condition  that 
they  could  generate  attenuation  events  of  at 
least  10  db.  In  this  way  not  only  a large 
portion  of  data  for  a given  day  were  neglected 
but  also  whole  days  were  not  considered,  so 
that  the  present  data  are  representative  of 
a period  of  time  much  larger  than  the  six 
days  over  which  they  were  collected.  Un- 
fortunately it  is  not  possible  to  establish 
exactly  the  representativeness  of  the  data 
because  not  enough  is  known  about  the  weather 
when  the  radar  wasn’t  in  operation  for  the 
reasons  above.  Further,  these  reasons  are 
not  independent  of  the  weather,  in  the  case 
of  power  failures,  thus  complicating  the 
problem  further.  In  this  work  it  has  been 
assumed  that  the  data  are  representative  of 
the  whole  sutmer,  although,  a distinct 
possibility  exists  that  significant  data  have 
been  missed. 

The  various  histograms  of  attenuation 
events  can  be  Interpreted  as  conditional  dis- 
tribution functions  with  the  condition  that 
the  probability  of  the  10  dB  event  is  given. 
They  were  integrated  to  give  probabilities 
and  normalized  in  such  a way  that  the  prob- 
ability of  the  10  dB  event  was  unity.  One 
example  is  given  in. fig.  2 for  a length  of 
5 miles  and  at  11.2  GHz. 

This  was  done  for  all  lengths  and 
all  frequencies.  To  obtain  absolute  prob- 
abilities it  is  enough  to  estimate  the  prob- 
ability of  the  10  dB  events  and  then  to 
scale  the  probabilities  for  higher  attenu- 
ation values  accordingly. 

As  exposed  in  chapter  3t  the  esti- 
mate of  the  probability  of  the  10  dBz  events 
can  be  obtained  from  raingauge  records,  thus 
allowing  construction  of  absolute  probability 
curves  for  all  attenuation  events.  The  ac- 
tual results  are  presented  in  paragraphs  2.5 
and  2.6. 


2.L  Radar  Calibration 

The  radar  equation  for  weather  radar 
with  STC  (sensitivity  time  control)  can  be 
written  as: 

where  Z is  the  reflectivity  factor,  is  the 
average  received  power  and  C is  the  calibra- 
tion constant  which  Includes  the  character- 
istics of  the  radar  and  the  dielectric  con- 
stant of  the  precipitation.  C changes  in 
time  and  has  been  determined  daily  following 
the  procedure  described  by  Smith  (9).  Tlie 
radar  data  have  then  been  used  to  calculate 
attenuation  statistics  for  5 mile  links  at 
11.2  GHz  and  these  have  been  compared  with 
the  one  obtained  by  a real  5 mile  link 
operating  at  Ormstown  within  radar  range  (10), 
fig.  2.  The  probability  curves  are  normalized 
to  1 for  10  dB  of  attenuation  and  are  obtained 
using  different  values  of  C in  steps  of  0.5  dB 
in  order  to  obtain  the  best  fit  between  the 
calculated  and  the  observed  statistics  (dotted 
line).  It  is  apparent  that  the  computed  sta- 
tistics are  very  sensitive  to  the  radar  cali- 
bration. The  probability  of  exceeding  30  dB 
changes  by  a factor  10  for  only  a 2 dB  vari- 
ation in  the  constant  C.  This  effect  is  in- 
fluenced by  hop  length  and  frequency  and 
would  be  different  if  the  curves  were  more 
logically  normalized  to  0 dB  but  this  example 
at  least  gives  an  idea  of  the  order  of  magni- 
tude. 2 dB  is  about  the  accuracy  one  can 
reasonably  expect  from  a careful  electrical 
calibration  of  the  radar.  The  importance  of 
a calibration  "ad  hoc"  is  then  apparent,  that 
is  a calibration  based  on  the  fit  of  the 
attenuation  statistics  obtained  from  a large 
sample  of  data.  In  this  case  the  discrepancy 
between  the  "electrical  calibration  constant" 
and  the  "attenuation  calibration  constant" 
was  about  U dB  which  compares  well  with  the 
value  of  3>J  found  by  1.  Strickland  in  Ottawa 
for  1970  data  (2),  following  a similar  method. 


2.5  Single  link  attenuation  probabilities 

Following  the  method  outlined  in  2.3 
and  2.4  attenuation  statistics  were  obtained 
from  radar  observations  in  the  summer  of  1970. 

Fig.  3 shows  the  probability  that  a 
given  attenuation  is  exceeded  at  11.2  GHz  for 
various  link  lengths.  For  comparison  fig.  3 
shows  the  probability  curve  obtained  from  the 
real  link  (dotted)  and  the  probability  curve 
obtained  for  5 mile  links  from  10  years  of 
raingauge  records  normalized  to  the  same  prob- 
ability for  10  dB  (dashed). 
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It  is  possible  to  observe  some  inter- 
esting features.  The  lowest  fraction  of  time 
observed  for  the  single  link  is  '*'10“°J  for 
10  link  years  is  '^10"7j  for  about  a total 
of  04  hours  of  radar  data  is  ‘^10*^.  All  the 
curves  are  exponential  until  a certain  value 
of  probability  and  then  they  follow  a law  of 
the  type  , The  single  link  is  exponen-  , 

tlal  until  -v3.10"5,  the  10  years  until 
and  the  radar  until  It  appears  tl.at 

increasing  the  sample  size  increases  pro- 
portionally the  range  of  probability  over 
which  the  process  observed  can  be  considered 
exponential.  This  suggests  that  the  process 
is  exponential  and  that  it  is  our  method  of 
observation  that  introduces  deviations  from 
the  exponential  form. 


Figs.  4 end  5 show  the  same,  statis- 
tics for  12.7  end  18.7  GHz.  The  hop  length 
dependence  is  shown  more  explicitly  in  figs. 
6t7»8»  where  the  probability  of  a given 
attenuation  is  plotted  versus  hop  length  at 
the  three  frequencies.  In  fig.  6 the  statis- 
tics from  10  years  of  raingauge  records  are 
plotted  for  comparison  (dotted  lines).  The 
fit  is  forced  on  the  10  dB  level. 


2.6  Joint  probabilities 

For  the  design  of  twin  route  radio 
systems  the  knowledge  of  Joint  probabilities 
of  parallel  links  is  required. 

The  knowledge  of  joint  probabilities 
of  adjacent  links  on  the  same  route  is  re- 
quired if  a too  conservative  design  is  to  be 
avoided  (11).  Both  these  probabilities  have 
been  evaluated  from  the  radar  data.  Flgs.9» 
1G,11,12,13|14|  show  joint  probabilities  for 
parallel  links  at  11.2  GHz,  for  values  of 
hop  length  of  2.3>  5i  7«5i  10|  12.5  and  15 
miles  versus  separation  distances.  Figs.l5» 
16,17»18,19»20  are  analogous  to  the  previous 
one  except  at  12.7*  Finally,  figs. 21,22, 23 
are  at  18.7  GHz  for  hop  length  of  2.5,  5> 

7«5  miles.  As  expected  the  joint  probabil- 
ities approach  the  probabilities  of  single 
links  when  the  separation  distance  approaches 
zero. 

Conditional  joint  probabilities  for 
links  on  the  same  route  are  given  for  ad- 
jacent links  and  links  separated  by  a link 
length.  It  is  preferred  to  use  conditional 
probabilities  Instead  of  absolute  joint 
probabilities  because  their  use  is  more 
direct  in  the  design  of  the  communication 
system. 

Figs.  24,25  show  these  statistics 
at  11.2  GHz,  figs.  26,27  at  12.7  and  figs. 
28,29  at  18.7  CHz. 


2.7  Accuracy  of  the  results 

There  are  many  causes  for  errors  in 
this  approach.  They  correspond  to  the  var- 
ious steps  necessary  to  arrive  at  the  final 
results  which  consist  of  an  estimation  of  the 
time  that  a given  attenuation  is  exceeded. 

For  this, measurements  of  reflectivity  are 
made.  The  errors  associated  with  these  mea- 
surements depend  on  various  factors,  some  of 
them  controllable,  some  not.  Random  phase 
fluctuations  are  a first  source  of  uncertainty 
which  can  be  reduced  increasing  the  number 
of  independent  samples  (Smith,  Marshall, 
Hitschfeld)  (4»6).  Reflectivity  gradients 
also  introduce  uncertainties  (Rogers)  (7)« 
Finally  the  reflectivity  is  stored  in  thres- 
holds, in  our  case  5 dB,  and  this  also 
introduces  errors  (Marshall)  (5).  Further, 
the  McGill  radar  makes  use  of  pulse  com- 
pression and  this  Introduces  other  uncertain- 
ties because  of  the  time-side-lobes  associated 
with  this  technique.  The  reflectivity  values 
stored  on  film  are  then  manually  transferred 
to  computer  cards,  thus  introducing  new  errors, 
although  the  careful  check  employed  limited 
this  to  a minimum.  The  simulation  of  attenu- 
ation introduces  more  errors  because  the 
exact  length  of  the  links  could  be  obtained 
only  in  particular  cases.  The  tolerance  on 
length  has  been  t 10^. 

The  effect  of  all  these  factors  is 
difficult  to  estimate.  For  this  reason  no 
clear  statement  on  the  accuracy  of  the  method 
is  available,  although  comparison  of  the  re- 
sults with  the  results  of  other  independent 
approaches  can  be  carried  on.  The  statistics 
obtained  with  this  approach  compare  very  well 
with  the  ones  obtained  by  a real  link  and  by 
the  raingauge  method,  figs.  3 and  6.  From 
these  comparisons  it  appears  that  the  accuracy 
of  the  method  is  at  least  satisfactory  be- 
cause it  produces  estimates  of  attenuation 
probabilities  that  agree  with  estimates  ob- 
tained independently. 


2.8  Conclusions 


Radar  provides  a powerful  tool  to 
evaluate  attenuation  statistics  for  microwave 
links.  The  main  advantage  over  other  methods 
is  the  enormous  sample  size  of  attenuation 
events  that  it  is  possible  to  collect  in  a 
few  days.  Of  the  six  days  considered  here, 
two,  27  July  and  10  August,  contributed  to 
most  of  the  events.  The  reason  is  that  one 
day  of  convective  activity  over  a region  of 
200  km  radius  offers  a full  range  of  vari- 
ability in  rainfall  rate  and  in  precipitation 
patterns. 

The  sample  size  depends  on  the 
sampling  time  of  the  data  and  on  the  areal 
coverage  of  the  radar.  With  the  present 
antenna  program  of  the  radar  the  minimum 
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i.a~pling  time  is  5 minutes.  It  was  observed 
on  a section  of  the  data  that  sampling  every 
IC  minutes  did  not  affect  the  results  so  that 
it  was  possible  to  reduce  the  data  analysis 
by  a factor  two.  It  was  also  found  that 
sampling  every  20  minutes  affects  the  re- 
sults. This  obviously  depends  on  the  area 
and  this  problem  will  be  studied  in  detail. 

Disadvantages  of  a radar  approach 
are,  among  others: 

- great  sensitivity  to  calibration,  which 
imposes  the  auxiliary  use  of  a microwave 
link  or  of  a ralngauge. 

- practical  impossibility  of  calculating, 
with  the  radar  alone,  absolute  prob- 
abilities. 

Furthermore,  year  to  year  variations 
in  the  statistics  obviously  cannot  be  extracted 
from  a few  days  of  radar  data,  but  they  can  be 
evaluated  by  ralngauge  records  or  directly  by 
microwave  link  data. 

The  joint  probabilities  presented 
here  contain  useful  information  about  the 
spatial  structure  of  precipitation  patterns. 

At  1-.7  GHz  rainfall  rate  and  specific  atten- 
uation are  almost  proportional,  so  that  a 
value  of  attenuation  is  proportional  to  the 
amount  of  rain  accumulated  along  the  same 
path  or  to  the  rain  accumulated  to  a point 
in  an  interval  of  time  equal  to  the  time  nec- 
essary for  t're  precipitation  pattern  to  travel 
a distance  equal  to  the  path  length,  “nit 
equivalence  of  space  and  time  statistics,  of 
precipitation  pattern  is  known  as  the  Taylor 
hypothesis.  Evidence  of  the  validity  of  this 
hypothesis  can  be  found  in  tlie  measurements 
of  autocorrelation  functions  cf  precipi ration 
patterns  by  Zawadzkl  (12)  and  in  the 
successful  application  of  this  hypothesis  to 
the  use  of  ralngauge  records  for  evaluating 
attenuation  statistics  (Chapter  3)»  From 
the  Joint  probabilities  of  attenuation  for 
parallel  or  adjacent  links  it  is  possible 
to  evaluate  the  joint  probabilities  of 
accumulated  precipitation  for  points 
separated  by  various  distances  and  for  the 
various  periods  of  accumulation.  This  sta- 
tistic is  of  interest  In  a variety  of  appli- 
cations . 


3* 

3.1  Introduction 

Rain  attenuation  statistics  can  be 
estimated  from  ralngauge  records  provided 
that  the  speed  of  translation  of  the  precip- 
itation patterns  is  known.  This  approach  is 
quite  appealing  because  suitable  ralngauge 
records  .ire  available  from  many  locations  in 
Canada  end  for  a sufficient  number  of  years. 
This  study  establishes  a technique  to  evaluate 
attenuation  statistics  from  the  ralngauge  re- 
cords. The  technique  Is  then  applied  to  cal- 
culate statistics  for  six  locations  across 
Canada  for  a period  of  time  of  ten  years  so 
that  year  to  year  variability  can  be  accounted 
for.  Further,  ten  years  of  data  provide  an 
adequate  sample  size.  In  this  chapter  the 
results  for  Montreal  are  presented,  the  re- 
sults for  the  ocher  locations  being  presented 
in  Chapter  4. 


3.2  Technique.  lescription  and 
testing 

The  method  of  employing  ralngauge 
records  to  geneiVite  attenuation  statistics 
is  based  on  the  concept  of  "synthetic  storm", 
first  used  in  1961  by  Hamilton  and  Marshall 
(13)  to  calculate  the  effects  of  rain  atten- 
uation on  data  obtained  by  a weather  radar 
operating  at  a wavelength  of  3 cm.  The  pre- 
sent method  is  an  application  of  the  Hamilton 
and  Marshall  concept  to  the  estiriatlon  of 
attenuation  over  a particular  ground  communi- 
cation link. 

As  a storm  or  rainfall  pattern  moves 
over  a ralngauge,  the  measured  rainfa'l  rate 
varies  with  time.  The  variations  arlsa  from 
two  effects,  the  advcction  of  the  spatial 
pattern  of  rain  and  changes  that  occur  within 
the  time  required  to  pass  over  the  ralngauge. 
These  effects  are  both  present  in  ralngauge 
records  and  are  not  separable  without  addition- 
al Information  about  the  structure  and  motion 
of  the  rain  pattern.  The  Haniilton-Marshall 
"synthetic  storm"  is  a description  of  the 
rain  pattern  in  terms  of  rainfall  rate  as  a 
function  of  distance  along  a line  in  the 
direction  of  storm  motion.  It  is  obtained 
from  a raingaugc  record  by  converting  the  raw 
data  of  rainfall  rate  versus  time  into  a 
function  of  distance  by  employing  a storm 
translation  velocity  to  transform  time  to 
distance.  In  the  original  Harallton-Marshal 1 
approach,  weather  radar  data  on  storm  motion 
provided  the  velocity  needed  for  the  tlme-to- 
distancf  conversion.  As  will  be  demonstrated 
bolc.w,  radar  data  may  not  be  essential  for 
establishing  this  velocity. 

Since  the  synthetic  storm  is  based 
or,  the  assumption  that  ail  time  cli.inges  ob- 
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served  at  a point  on  the  ground  arise  from 
the  advectlon  of  spatial  variations,  It  will 
not  give  an  exact  description  of  the  distri- 
bution of  rain  with  distance.  It  may  be, 
however,  that  statistical  properties  of  a 
large  number  of  synthetic  storms  are  approxi- 
mately the  same  as  the  corresponding  statis- 
tical properties  of  real  storms.  This  is  the 
hypothesis  tested  In  the  present  study,  where 
the  properties  of  Interest  are  certain  Inte- 
grals of  the  rainfall  rate  over  specified 
distances. 

This  hypothesis  is  also  known  as 
Taylor’s  hypothesis.  Evidence  of  It's  valid- 
ity can  be  found  not  only  In  this  work  but 
also  in  Zawadzki  (12). 

The  synthesized  storms  express  rain- 
fall rate  R as  a function  of  distance  x. 

As  the  first  step  in  this  analysis,  rainfall 
rate  is  converted  to  specific  attenuation  Y 
through  the  formula 

Y •=  kR**  (dB/mi)  (3 •2.1) 


where  k and  <*•  are  empirical  parameters  depend- 
ing on  radio  frequency.  The  attenuation  over 
a distance  L ma^  be  calculated  by  Integrating 
Y(x)  over  an  interval  of  length  L: 

A(x^,L)  - I Y(x)dx  . (3.2.2) 

X 

o 

The  result  of  this  integration  will  generally 
depend  or.  the  starting  point  Xg  as  well  as 
the  length  of  the  Interval.  The  variation 
of  attenuation  with  time  can  then  be  simu- 
lated by  changing  Xg  at  a rate  equal  to  the 
velocity  of  storm  notion.  In  order  to  test 
the  technique  the  measurements  obtained 
from  the  BNR  link  In  Ormstown  (fig.  l)  were 
used.  Three  ralngauges  were  located  along 
the  link  and  they  supplied  the  rainfall  data. 

Fig.  30  Illustrates  a particular 
synthetic  storm  in  terms  of  both  rainfall 
rate  and  specific  attenuation,  with  the  time- 
to-distance  conversion  Indicated. 

Attenuation  was  calculated  as  a 
function  of  time  for  this  storm  by  use  of 
(3.2.2),  with  L set  at  5 nil  and  Xg  changing 
at  a rate  of  14  ml/hr  (as  obtained  by  radar). 
Values  of  k and  ?(.  were  selected  to  corres- 
pond to  11.2  GHz.  The  result  of  this  calcu- 
lation is  shown  In  Fig.  31  and  compared  with 
the  actual  attenuation  observed  on  the  link. 

That  the  calculated  and  actual  atten- 
uations are  quite  different  Is  not  surprising. 
First,  the  calculated  attenuation  cannot  ac- 


count for  any  time  variations  occurring  In 
the  storm  as  it  moves  along.  Secondly,  the 
simulated  attenuation  record  corresponds  to 
a link  that  Is  parallel  to  the  direction  of 
storm  motion.  In  this  example  (as  In  most 
of  the  actual  cases)  the  storm  had  an  appre- 
ciable component  of  motion  across  the  Omstown 
link.  Fig.  31  also  shows  that  the  simulated 
attenuation  has  a smoother  curve  than  the  ob- 
sen^ed  attenuation.  This  is  associated  with 
the  fact  that  the  ralngauge  records  do  not 
allow  rainfall  rate  to  be  measured  with  a 
time  resolution  better  than  about  1 min.  No 
attempt  was  made  to  calculate  attenuation 
with  a finer  resolution. 

The  experiment  has  been  carried  on 
for  the  summers  of  1971  and  1972  and  statis- 
tics from  the  link  and  the  gauges  are  com- 
pared in  fig.  32. 

In  these  comparisons  a set  of  values 
of  k and  m has  been  chosen  so  that  the 
best  fit  is  obtained.  The  values  of  k and  cx 
are  given  In  table  1 for  the  three  frequencies 
of  Interest,  and  they  compare  well  with  the 
values  given  in  the  literature,  (14)  • 

In  spite  of  the  fact  that  an  individual  syn- 
thetic storm  may  not  lead  to  an  attenuation 
record  that  agrees  with  observations,  fig. 

32  Indicates  a good  agreement  provided  that 
the  sample  size  Is  sufficient.  While  It 
may  have  been  expected  that  the  method  of 
calculating  attenuation  would  lead  to  some 
systematic  bias  effect,  since  all  the  simu- 
lated links  are  in  the  direction  parallel  to 
storm  motion,  this  does  not  appear  to  be  the 
case.  Moreover,  the  results  here  indicate 
that  the  1 minute  resolution  of  the  ralngauges 
is  adequate,  at  least  over  paths  as  short  as 
5 miles. 


Dependence  on  the  parameters 


Velocity.  The  translation  velocity 
of  the  storm  plays  an  essential  role  In  this 
analysis.  Radar  provides  the  best  estimate 
of  this  velocity.  But  because  radar  facili- 
ties are  not  as  common  as  ralngauges,  an 
attempt  has  been  made  to  estimate  the  trans- 
lation velocity  from  conventional  radiosonde 
data.  It  has  been  found  that  the  700  mb  wind 
nearest  in  time  is  a good  estimate  of  the  ve- 
locity. 700  mb  corresponds  on  average  to 
3 km  above  sea  level  and  is  often  referred 
to  as  "steering  level".  Fig.  33  compares 
probabilities  obtained  using  radar  and  TOO  rob 
velocities.  The  soundings  are  from  Manlwakl 
120  miles  away  from  the  experimental  location 
(fig.  1).  This  comparison  is  considered 
successful  and  it  Justifies  the  use  of  radio- 
sonde velocities  as  estimates  of  the  trans- 
lation velocity  of  precipitation  patterns  - 
k and  . 
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The  probability  of  a given  value  of 
attenuation  depends  on  the  value  of  k and  « . 
In  this  study  no  direct  evaluations  of  k and 
verc  available  so  that  values  were  chosen 
In  order  to  have  the  best  fit  between  calcu- 
lated and  observed  statistics.  This  Is  a 
circular  argument.  To  increase  the  confidence 
in  Che  results  the  sensitivity  of  the  results 
on  k and  > has  been  determined.  The  prob- 
ability Chat  a given  attenuation  is  exceeded 
car.  be  considered  exponential: 


P(a)  P e 
o 


“A/  Ar 


(3.3.1) 


where  ^isa  function  of  hop  length  and  Ao  is  a 
function  of  frequency  and  hop  length;  A is 
dependent  on  k and  ».  . 

AP 

It  is  of  interest  In  function  of 
Ak  and  A o>.  . 

% doesn’t  depend  on  k and  <>■  because 
it  is  Che  percentage  of  time  when  it  rains. 
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3. A Ten  years  statistics  - single 
links  


Using  Che  technique  described  In  3.2, 
ten  years  of  ralngauge  records  for  Montreal 
have  been  analyzed.  The  years  go  from  1961 
to  1970  and  each  year  goes  from  May  to 
September. 

Fig.  3A  shows  Che  probabilities  ob- 
tained for  5 miles  hop  length  at  11.2  GHz 
for  the  10  different  years.  196A  has  a 
peculiar  behavior:  one  storm  is  such  that 

it  extends  the  curv’e  up  to  7.-  dB.  This  is 


the  record  storm  for  Montreal  with  an  Inten- 
sity of  A60  nm/hr.  The  return  period  of  this 
event  exceeds  100  years  so  chat  it  cannot  be 
Included  In  averages  over  a ten  year  period. 
For  this  reason  it  has  been  disregarded  In 
Che  analysis. 

It  is  not  clear  Co  the  author  which 
Is  the  best  form  of  presenting  statistics  of 
events  like  the  one  In  fig.  34.  It  has  been 
decided  to  plot  average  probabilities  over 
ten  years  for  various  hop  lengths  and  frequen- 
cies (in  this  paragraph)  and  to  describe  year 
Co  year  variability  separately  (3.8).  An 
alternative  approach  is  Illustrated  in  fig. 

35,  where  not  only  average  probabilities  are 
plotted  but  also  probabilities  corresponding 
Co  certain  percentages  of  the  population. 

The  probabilities  of  a given  attenuation  In 
Che  various  years  are  assumed  to  be  random 
events  normally  distributed.  Although  this 
approach  Is  limited  here  Co  one  frequency 
and  one  link  length  It  can  be  extended  to 
other  values  using  the  data  on  year  Co  year 
variability  in  3«8. 

Figs.  36,37,38  plot  probabilities 
versus  attenuation  at  11.2,  12.7,  18.7  GHz; 
for  hop  lengths  of  2.5,  5,  7.5,  10,  12.5, 

15  miles.  The  curves  are  straight  lines 
until  a certain  percentage  of  time  ('^3.10"° 
of  one  year)  so  that  Che  probabilities  can 
be  considered  exponential  until  that  prob- 
ability level.  Comparison  with  radar  data 
for  larger  sample  size  (paragraph  2.5) 
suggests  that  this  Is  an  observation  effect 
due  to  insufficient  sample,  so  Chat  Che 
probabilities  should  be  considered  exponen- 
tial all  Che  way. 


3.5  Joint  probabilities 

With  the  same  technique  described 
in  3.2  it  it  possible  Co  calculate  statistics 
of  Joint  attenuation  events  for  successive 
links,  these  statistics  being  of  Interest  In 
the  planning  of  radio  systems.  Statistics 
for  adjacent  links  are  plotted  in  figs.  39, 
40,41  In  the  form  of  conditional  probabil- 
ities, the  condition  being  that  a given 
attenuation  is  exceeded  on  one  of  the  links. 

At  18. 7 OHz  Che  plots  are  limited  to  a length 
of  5 miles.  Figs.  42,43,44  present  the  same 
statistics  for  links  separated  by  a hop  length. 


3.6  Durations 


Table  II  presents  a comparison  of 
the  durations  of  Che  attenuation  events  In 
the  1971  Ormstown  experiment.  It  appears 
that  Che  link  and  Che  ralngauges  agree  al- 
though Che  sample  size  Is  so  small  Chat  It 
is  impossible  to  assess  the  culpability  of  the 
ralngauge  to  reproduce  the  statistics  of 
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durations  experienced  by  the  link.  Neverthe- 
less | it  appears  that  the  two  sets  of  obser- 
vations are  consistent,  i'lstograms  of  dura- 
tions of  fade  are  presented  for  various 
lengths  at  11.2  CHz  only  in  the  tables  III, 
IV,  V,  VI,  VII,  VIII.  As  an  example,  figs . 

45  and  L.i  show  the  dependence  on  hop  length 
and  on  attenuation  level. 


3.7  Frepuency  dependence 

From  the  10  years  statistics  presented 
in  figs.  36,37,38  it  is  possible  to  obtain 
some  information  about  frequency  dependence 
of  attenuation  probabilities.  In  figs.  47,48, 
49  attenuation  of  events  with  the  same  prob- 
ability and  hop  length  but  different  frequency 
are  correlated  so  that  the  following  empirical 
results  are  found: 


18.7  " ^11.2 

(3.7.1) 

18.”  " ^12.7 

(3.7.2) 

12.-  ^*25  ^11.2 

(3.7.3) 

which  relate  the  values  of  attenuation  at  two 
frequencies  for  the  same  length  and  probability. 

From  the  above  relations  one  can  ob- 
tain: 

* (f)'-’'’  0.7.0 

*^fo  0 

which  at  least  in  the  range  11  to  IS  GHz  allows 
one  to  calculate  empirically  the  attenuations  of 
an  event  at  frequency  f given  the  attenuation 
of  the  sane  event  at  f^,,  the  same  event  meaning 
the  same  link  length  and  probability. 


3.3  Year  to  year  variability 

As  can  be  seen  in  fig.  34  there  is 
a substantial  variation  in  the  probability 
of  a given  attenuation  event  across  the  years. 

For  any  length  and  frequency  the 
probability  that  a given  attenuation  is  ex- 
ceeded for  a given  year  is  very  likely  to  be 
a Gaussian  variate.  tests  for  typical 
values  cannot  reject  this  hypothesis  at  any 
probability  level  but  the  number  of  samples 
is  so  small  that  the  significance  of  the  test 
is,  at  least,  dubious. 

These  yc.irly  variations  have  to  be 
taken  into  account  in  the  planning  of  radio 
systems  and  also  in  establishing  the  duration 
of  experiments  to  collect  attenuation  statis- 
tics. These  variations  depend  on  frequency 
and  i ink  length  in  a com.pl  ic.ited  vay.  In 


fig.  50  the  coefficient  of  variation  of  the 
probability  of  a given  event  is  plotted  as 
a function  of  hop  length  for  discrete  atten- 
uation levels.  The  coefficient  of  variation 
is  defined  as  follows: 

g 

c.o.v.  — — X 100 


where  S is  the  sample  variance  and  ^ is  the 
sample  mean,  the  variable  being  the  prob- 
ability of  a given  event. 


4.  ATTENUATION  BY  RAINGAUGE  ACROSS  CANADA 
4»1  Introduction 

For  the  purpose  of  obtaining  some 
information  about  the  geographical  dependence 
of  the  attenuation  statistics  of  interest,  5 
more  locations  have  been  selected  to  carry 
on  a study  similar  to  the  one  carried  on  for 
Montreal.  Fig.  51  shows  these  locations  and 
the  radiosonde  stations  whose  upper  air  data 
were  employed  in  the  analysis.  As  in  Montreal, 
10  years  of  raingauge  data  were  analyzed  for 
each  location  covering  the  months  from  M.ay  to 
September  for  the  decade  from  1962  -o  1971. 

The  actual  locations  were  chosen  as  corres- 
ponding to  highly  populated  regions  and  with 
the  intention  of  covering  the  possible  routes 
of  Canadian  radio  systems.  The  west  coast 
and  the  Rocky  Mountains  were  excluded  be- 
cause of  the  strong  '•i.ographic  character  of 
the  precipitation  there. 

The  results  of  this  study  are  cer- 
tainly not  conclusive  or  complete,  leaving 
large  parts  of  the  country  uncovered.  In 
particular,  the  region  between  Halifax  and 
Montreal  needs  further  analysis  as  well  as 
the  region  between  Edmonton  and  Winnipeg. 
Nevertheless  an  overall  pattern  appears  from 
the  available  data.  There  is  a region  con- 
taining Montreal,  Toronto,  Thunder  Bay  and 
Winnipeg  where  attenuation  statistics  are 
severe  and  particularly  similar  for  the  four 
locations.  East  of  this  region,  maritime 
climate  generates  less  severe  statistics 
(Halifax),  and  westward  (Edmonton)  less 
severe  statistics  occur  in  correspondence  to 
a continental  climate.  The  boundaries  of 
this  region  are  not  well  defined  at  this 
stage  of  the  analysis. 

In  carrying  on  this  study  a further 
hypothesis  has  to  be  introduced.  That  is 
that  the  700  mb  winds  are  representative  of 
the  precipitation  patterns  motion  all  over 
the  country.  This  is  largely  arbitrary  be- 
cau«^c  this  hypothesis  b.is  been  tested  only 
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in  Montreal  (flg«  33 )•  Nevertheless  this 
appears  to  be  a reasonable  hypothesis  and 
certainly  not  a critical  one}  fig.  52  shows 
Winnipeg  statistics  calculated  using  radio- 
sonde data  from  International  Falls  and 
Sault  Ste.  Marie.  The  curves  are  very  close 
to  each  other  thus  Insuring  a certain  re- 
liability to  the  study.  In  the  following 
paragraphs  the  statistics  for  the  various 
locations  are  presented. 


4.2  Halifax 

Figs.  53f54i55  show  proiabllltles 
versus  attenuation  at  11. 2,  12. 7>  18.7  GHz. 

Like  for  Montreal  the  probability  :hat  a 
given  attenuation  Is  exceeded  decreases 
exponentially  versus  attenuation  at  least 
to  a probability  of  about  3>10~^  of  one  year. 

The  Joint  probabilities  (figs.  56t 
57»58,59t60|61)  are  very  small  for  every  value 
of  attenuation,  length  and  frequency.  This 
result  has  Important  consequences  for  the 
design  of  radio  systems  because  attenuation 
events  on  consecutive  links  can  then  be  con- 
sidered Independently. 

Tables  IX  to  XIV  present  histograms 
of  duration  at  11.2  GHz  for  various  hop  lengths. 
Fig.  62  shows  the  coefficient  of  variation 
across  ten  years.  For  the  definition  of  these 
curves  sec  paragraph  3*8.  The  cun’es  in  fig. 

62  di'ft-  from  the  Montreal  one.  fig.  50,  in- 
the  sense  that  the  variability  In  Halifax  Is 
substantially  larger,  and  the  length  depen- 
dence Is  different,  suggesting  smaller  high 
rate  rain  cells  than  In  Montreal. 


4.3  Toronto 

Figs.  63»6~,65  show  probabilities 
versus  attenuation  at  11.2,  12.7,  18.7  GHz. 
The  curves  are  very  close  to  those  for 
Montreal  although,  as  will  be  seen  In  para- 
graph 4.7,  the  precipitation  statistics  arc 
different. 

The  joint  probabilities  (figs.  66, 
67i68,69,70,7l)  suggest  complicated  Inter- 
actions between  cell  size  and  hop  lengths. 
They  are  generally  higher  than  In  Halifax 
and  similar  to  those  in  Montreal.  The  values 
are  such  to  make  these  statistics  an  Impor- 
tant factor  In  the  design  of  radio  systems. 

Histograms  of  duration  at  11.2  GHz 
arc  given  In  the  tables  XV  to  XX  in  corres- 
pondence to  various  lengths. 

The  coefficients  of  variations  are 
plotted  in  fig.  72  (for  the  definition  sec 
paragraph  3*8)*  It  appears  that  the  yearly 
variability  is  generally  somehow  higher  than 
in  Motitroil  and  less  than  in  Halifax. 


4.4  Thunder  Bay 

As  the  name  says.  Thunder  Bay  Is 
quite  stormy  at  least  for  what  concerns  atten- 
uation. Figs.  73*74,75  show  probabilities 
versus  attenuation  at  11.2,  12.7*  16«7  GHz. 
These  curves  are  very  close  to  the  ones  for 
Montreal  and  Toronto.  The  joint  probabilities 
(figs.  76*77*78*79*80,81)  are  close  to  the 
Montreal-Toronto  ones  for  short  lengths  and 
somehow  less  for  longer  lengths  suggesting 
smaller  sizes  of  rain  cells.  Also  In  Thunder 
Bay  the  values  of  joint  probabilities  are 
such  to  make  them  relevant  in  the  design  of 
radio  systems. 

Histograms  of  duration  are  presented 
In  Tables  XXI  to  XXVI.  The  coefficients  of 
variation  are  plotted  in  fig.  82  (for  the 
definition  see  paragraph  3*8).  Relative  to 
the  other  locations  the  variability  In 
Thunder  Bay  Is  higher  than  in  Montreal, 
slightly  higher  than  in  Toronto  and  lower 
than  In  Halifax. 


4.5  Winnipeg 

Figs.  83,84*85  show  probabilities 
versus  attenuation  at  11.2,  12.7,  18.?  GHz. 
Again  these  curves  are  very  similar  to 
those  for  Montreal.  In  calculating  these 
curves  some  extreme  events  were  discarded 
from  the  data. 

The  joint  probabilities  (figs.  86, 
87,88,39,90,91)  are  very  similar  to  those 
for  Montreal. 

Histograms  of  durations  are  presented 
'in  Tables  XXVII  co  XXXII.  The  coefficients 
of  variation  are  plotted  in  fig.  92  (for  the 
definition  see  paragraph  3*8).  Relative  to 
other  locations  the  variability  in  Winnipeg 
Is  higher  except  for  Halifax. 


4.6  Edmonton 

A well  established  pattern  of  simi- 
larities in  the  attenuation  statistics  Is 
encountered  moving  westward  from  Montreal. 

This  pattern  presents  a marked  discontinuity 
In  Edmonton  and  whether  this  is  the  result 
of  a gradual  change  or  of  an  abrupt  discon- 
tinuity somewhere  between  Winnipeg  and  Edmonton 
Is  not  yet  known.  The  available  climato- 
logical data  do  not  furnish  an  explanation 
of  this  problem  and  It  will  probably  be 
necessary  to  analyze  ralngauge  data  for  lo- 
cations between  Winnipeg  and  Edmonton. 

Figs.  93,94  and  95  show  probabil- 
ities versus  attenuation  at  11.2,  12.7  and 
18.7  GHz.  For  any  length  and  frequency  the 
probability  of  exceeding  a given  attenuation 
is  much  le.s:.  than  in  any  other  locality  in 
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this  study,  thus  allowing  much  longer  hop 
lengths. 

The  Joint  probabilities  (figs.  96, 
97,^8,99,100  and  101}  are  so  small  that  no 
advantage  can  be  obtained  by  taking  them  In- 
to account  for  the  design  of  radio  systems. 

The  behaviour  of  the  single  link  and  joint 
link  statistics  in  Edmonton  can  be  interpreted 
in  terms  of  frequency  of  occurrence  of  given 
rain  events  and  their  size.  Something  in 
this  line  will  be  attempted  in  the  next 
paragraph. 

Histograms  of  durations  at  11.2  GHz 
are  given  in  Tables  XXXIII  to  XXXVIII.  The 
coefficients  of  variation  are  plotted  In 
fig.  102  (for  the  definition  see  paragraph 
3.8),  they  are  generally  higher  than  in  the 
other  locations  but,  of  course,  the  smaller 
sample  .-:ize  of  the  events  has  to  be  taken 
into  account  in  explaining  this  behaviour. 


Comp^arlsons 


As  stated  in  para.  4>1,  the  six 
locations  considered  can  be,  in  many  respects, 
distributed  in  three  groups  with  different 
properties.  Fig.  103  shows,  the  attenuation 
exceeded  for  a probability  of  10'^  of  one 
year  in  function  of  hop  length  for  the  six 
locations.  It  is  apparent  that  Montreal, 
Toronto,  Thunder  Bay,  Winnipeg  behave  very 
much  In  the  same  way,  while  Halifax  and 
Edmonton  have  a different  behaviour.  A 
similar  pattern  could  be  observed  for  the 
joint  probabilities  although  less  clearly  be- 
cause these  quantities  are  determined  with 
a smaller  sample  size  and  thus  with  less 
statistical  accuracy.  The  year  to  year 
variability  doesn't  follow  a pattern  as 
clear  as  the  one  for  the  probabilities,  and 
a conparlson  between  the  locations  must  take 
into  account  the  relative  sam.ple  sizes. 

For  these  reasons  some  statistics  of  the 
precipitation  processes  In  the  various  loca- 
tions arc  presented  here,  having  been  ob- 
tained from  the  same  raingauge  records  used 
for  the  attenuation  studies. 

Table  XXXIX  shows  the  number  of 
precipitation  events  in  ten  years  exceeding 
a given  rainfall  rate  in  the  six  locations, 
and  Table  XL  shows  the  features  of  the 
population  of  the  event  sizes.  It  appears, 
for  instance,  that  the  behavior  of  Edmonton 
statistics  is  consistent  with  these  data 
which  show  rarer  and  shorter  events  which 
corresponds  to  smaller  probability  of  a 
given  attenuation. 


/,.8  Planning  a radio  system  across 
Canada 

The  result  presented  here  so  far 
can  be  used  for  the  planning  of  a radio 
system  across  Canada,  at  least  in  a general 
sense.  They  are  by  no  means  sufficient  for 
a practical  design  for  which  a more  detailed 
coverage  of  the  country  would  be  needed.  In 
a general  sense  the  data  presented  here  allow 
some  considerations: 

It  appears  that  there  is  one  large 
region,  at  least  extending  from  Montreal  to 
Winnipeg,  where  the  probability  of  exceeding 
a given  attenuation  is  almost  uniform  so 
that  a uniform  criterion  of  design  is  allowed. 
Inside  this  region  year  to  year  variability 
is  larger  in  certain  places  than  in  other;i 
so  that  in  order  to  have  uniform  performances 
some  allowances  have  to  be  made  in  this  res- 
pect. In  this  region  the  spatial  structure 
of  the  precipitation  pattern  is  such  that 
joint  events  of  sizeable  attenuations  are 
frequent  enough  to  make  it  worth  taking  them 
into  account  for  more  economical  planning. 

Outside  this  region,  both  eastward 
and  westward,  the  probabilities  are  lower 
thus  making  possible  the  use  of  longer  hop 
lengths  although  this  advantage  is  partially 
compensated  by  larger  year  to  year  variabil- 
ity and  independency  of  events.  The  boun- 
daries of  the  region  are  not  known,  so  that 
it  is  Impossible  to  evaluate  accurately  the 
cost  of  a system.  For  this  further  studies 
arc  necessary. 
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Fig.  1.  The  experimental  link  was  located  In  Ormstown 
20  miles  south  of  Montreal. 


12< 


RELATIVE  PROBABILITY 
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10  20  30  40  A (dB) 

Fig.  2.  Relative  probabilities  that  a given  attenuation 
is  exceeded  (continuous  lines)  as  obtained  for  different 
values  of  the  calibration  constant  (in  a relative  scale). 
The  dotted  curve  is  obtained  from  the  experimental  link. 
All  curves  are  normalized  to  1 for  10  dB  of  attenuation. 
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Figs.  6»7»8  (Left;  Below  Left;  Below  Right) 
Probability  that  a given  attenuation  is 
exceeded  at  11. 2»  12. 7»  end  18.7  CHz  ver- 
sus hop  length.  In  Fig.  6 the  dotted 
lines  are  from  10  years  of  raingauge  re- 
cords. A fit  has  been  forced  at  the 
10  dB  level.  The  scale  on  the  left  is 
fraction  of  one  year. 
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Figs.  21,22,23  (Left;  Below  Leftj  Below 
Right).  Probability  that  a given  atten 
nation  is  Jointly  exceeded  on  parallel 
links  at  18. 7 GHz  versus  separation  dls 
1^-4  tance  for  various  hop  lengths. 
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Fig.  24*  Conditional  probability 
that  a given  attenuation  is  ex- 
ceeded Jointly  on  consecutive 
links  versus  hop  length  at  11 >2 
GHz. 
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Fig.  25*  Conditional  probability 
that  a given  attenuation  is  ex- 
ceeded Jointly  on  links  separated 
by  a hop  length  ver,.u8  hop  length 
at  11.2  GHz. 
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Fig.  26.  Conditional  probability 
that  a given  attenuation  is  ex- 
ceeded Jointly  on  consecutive 
links  versus  hop  length  at  12.? 
GHz. 
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Fig.  27.  Conditional  probability 
that  a given  attenuation  is  ex- 
ceeded Jointly  on  links  separated 
by  a hop  length  versus  hop  length 
at  12.7  GHz. 


Fig.  2S.  Conditional  probability  that  a 
given  attenuation  ia  exceeded  Jointly 
on  consecutive  links  versus  hop  lengths 
at  18.7  GHz. 
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Fig.  29.  Conditional  probability  that  a 
given  attenuation  is  exceeded  jointly 
on  links  separated  by  a hop  length 
versus  hop  length  at  18.7  GHz. 
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Fig.  30  (Above).  Synthetic  storm  of  the  20th  August  1971  (continuous  line).  The 
dotted  line  refers  to  specific  attenuation  at  11.2  GHz. 

Fig.  31  (Below).  The  storm  of  Fig.  30  has  been  used  to  calculate  the  attenuation 
for  a 5-mile  link  at  11.2  GHz.  For  comparison  the  actual  attenuation  experienced 
by  this  link  is  plotted.  Note  the  coincidence  of  the  peak  at  22  minutes.  This 
is  due  to  chance  because  the  storm  hit  the  raingauge  with  its  maximum  intensity. 
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Fig.  32.  Probability  that  a 
given  attenuation  is  exceeded 
for  a 5-mile  link  at  11.2  GHz 
as  calculated  from  raingauge 
(dashed)  and  measured  on  the 
real  link  (continuous).  The 
data  are  for  I97I  and  1972. 


Fig.  33*  1971  statistics  using 

different  estimates  of  the 
velocity  of  precipitation 
patterns.  The  curve  with 
two  dots  is  obtained  using 
one  value  of  velocity  for 
the  entire  summer.  This  is 
to  show  the  relative  insensi- 
tivity of  the  results  to  the 
velocity  used  in  the  calcu- 
lations . 
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Fig>  34*  Probabilities  that 
a given  attenuation  is  ex- 
ceeded at  11.2  GHz  for  5~ 
mile  link  for  10  years  in 
Montreal. 


fig*  35*  Probability  that  a 
given  attenuation  is  ex- 
ceeded at  11.2  GHz  for  5> 
mile  link  for  a lO  year 
period.  The  curves  labelled 
0.9  and  C.9999  correspond 
to  \Q%  and  .01^  of  the  popu- 
lation assumed  to  be  gaussian. 
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Figs*  36,37,38  (Top  Leftj  Bottom  Left;  Above).  Probability  that  a given 
attenuation  is  exceeded  at  11.2,  12.7,  18.7  GHz  for  various  hop  lengths. 
The  curves  are  obtained  averaging  10  years  of  data. 
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Figs.  39»40,41  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  Is  exceeded  contem- 
poraneously on  adjacent  links  at  11.2, 
12.7,  and  18.7  GHz. 
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Fig.  45*  Histograms  of 
durations  of  25  dB 
i.ismiits  attenuation  events  at 

11.2  GHz  for  various 

j I 1 — hop  lengths  in  10  years 
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Fig*  46.  Histograms  of 

@ ti  2 GHi  durations  of  events  at 

L'Smi  a given  attenuation 

5 for  a 5 mile  link  at 

11.2  GHz  in  10  years. 
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Figs.  47»48,49  (LeftJ  Below  Left;  Below 
Right).  Correlogram  of  attenuation 
events  at  11.2,  12. 7»  18.7  GHz.  The 
events  have  the  same  probability  level 
and  are  for  the  same  hop  length. 
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Fig.  50.  Coefficient  of  variation  of  th^ 
is  exceeded  versus  hop  length  for  11.2 


Fig.  51  • Map  of  Canada  showing  the  location  of  the  raingauge  stations  and  rawlnsonde 
stations  considered  in  this  study. 


The  two  curves  are  calculated  using  wind 
velocities  from  International  Falls  (con- 
tinuous) and  Sault  Ste.  Marie  (dotted). 
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Figs.  53f54»55  (LeftJ  Below  Left;  Below 
Right).  Probability  that  a given 
attenuation  is  exceeded  at  11. 2»  12.7 
and  18.7  CHz  for  various  hop  lengths. 
The  curves  are  obtained  averaging  10 
years  of  data. 


COND-  PROB.  FOR  ADJACENT  LINKS  COND  PROD.  FOR  ADJACENT  LINKS 


HOP  LENGTH 


Figs.  56f57»58  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  is  exceeded  contem- 
poraneously on  adjacent  links  at  11. 
12.7  and  18.7  GHz. 


C 2 8 5 7 8 10  12  8 IS  mile*  012345  mil»» 

HOP  length  hop  length 


r-.f" 

».'0 


< 


COmO  HrtOJ.  FOI(  LINKS  SEPARATED  BT  A HOP  LENSTH  ^ CONO  PROf)  FOR  LINKS  SEPARATED  BY  A MOP  LENGTH 


T 

t (D  12  r GMi 


malifak 


15  mil«« 


coNo  PROD,  ron  links  sePARATto  er  a hop  lencth 


Figs.  59*60,61  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  is  exceeded  contem- 
poraneously on  two  links  separated  by  a 
link  length  at  11.2,  12.7  and  18.7  GHz. 
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Fig.  62.  Coefficient  of  variation  of  the  yearly  probability  that  a given  attenuation 
is  exceeded  versus  hop  length  for  11.2,  12*7  and  18.7  GHz. 
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Figs.  66»67«68  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  Is  exceeded  contem- 
poraneously on  adjacent  links  at  11.2, 
12.7  and  18.7  GHz. 
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Fig.  72.  Coefficient  of  variation  of  the  yearly  probability  that  a given  attenu 
ation  is  exceeded  versus  hop  length  for  11. 2»  12.7  end  18.7  GHz. 
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Figs.  73»74,75  (Left;  Below  Left;  Below 
Right).  Probability  of  exceeding  a given 
attenuation  at  11.2,  12.7  and  18.7  GHz 
for  various  hop  lengths.  The  curves  are 
obtained  averaging  10  years  of  data* 
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Figs.  76,77*78  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  is  exceeded  contem- 
poraneously on  adjacent  links  at  11.2, 
12.7,  and  18.7  GHz. 
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Flg«.  79,80,81  (L«ftj  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  la  exceeded  contem* 
porancously  on  two  llnka  aeparatcd  by  a 
hop  length  at  11.2,  12»7  ®nd  18.7  GH*.  ^ 
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COMO.  PnOB.  FOR  ADJACENT  LINKS  COKD  PilOB.  FOR  ADJACENT  LINKS 


Figs.  86)87*88  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  la  exceeded  contem- 
poraneously on  adjacent  links  at  11.2, 
12.7  and  18.7  GHz. 
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Pigs.  89»90»91  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  is  exceeded  contem- 
poraneously on  two  links  separated  by 
a hop  length  at  11.2,  12.?  and  18.? 

GHz. 
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Fig.  92.  Coefficient  of  variation  of  the  yearly  probability  that  a given  attenuation 
is  exceeded  versus  hop  length  for  11. 2«  12.7  and  18.7  GHz. 
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Flg».  96,97»98  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
given  attenuation  is  exceeded  contem- 
poraneously on  adjacent  links  at  11.2, 
12.7  end  18.7  GHz. 
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Figs.  99>100»101  (Left;  Below  Left;  Below 
Right).  Conditional  probability  that  a 
® " 2 GH<  EDMONTON  glven  attcnuatlon  is  exceeded  contem> 

poraneously  on  two  links  separated  by  a 
link  length  at  11. 2»  12.7  and  18.7  GHz. 
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Fig.  103*  Attenuation  exceeded  at  18.7  GHz 
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function  of  hop  length  for  the  various 
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TABLE  mix 


Number  of 

rain  events 

in  10  years 

for  various 

rainfall 

* 

rates 

R (mm/hr ) 

10 

30 

50 

70 

90 

Halifax 

483 

159 

66 

39 

26 

Montreal 

699 

335 

123 

76 

54 

Toronto 

423 

188 

105 

70 

50 

Thunder  Bay 

436 

210 

99 

68 

53 

Winnipeg 

460 

179 

92 

60 

41 

Edmonton 

286 

94 

38 

14 

11 

An  event  of  a given  rainfall  rate  happens  when  rainfall  rate 
exceeds  the  given  value.  Unfortunately,  in  the  same  precipitation 
pattern  this  can  happen  more  than  one  time. 


TABLE  XL 


Statistical  properties  of  the  sizes  of  rain  events  for 
various  rainfall  rates.  For  each  rainfall  rate  three 
lengths  are  given  in  km.  They  correspond  to  the  30^ 
level  (median  size),  the  70^  level  and  the  90%  level.* 
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<.5 

<.5 
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13.5 
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1.5 
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2.1 

1.3 

1.7 

1 

.7 

3.5 

2.5 

3.5 

2.3 

2 

8.5 

6.5 

7.3 

5.5 

3.7 

Toronto 

2.75 

1.7 

1.1 

1 

? 

4.8 

3.4 

2 

1.5 

1.3 

11.5 

7.2 

5.5 

3.5 

3.3 

Thund«;r  Bay 

2.5 

1.5 

1.2 

1 

.75 

5 

3.5 

2.2 

1.4 

1.1 

12.5 

7.5 

6.5 

5 

2.8 

Winnipeg 

2.4 

1.5 

1.2 

.75 

.7 

4 

2.5 

2 

1.25 

1 

10 

5.5 

4.5 

2;  5 

2.5 
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2.2 

1.2 

1 

0.5 

<.5 

3.9 

2.6 

1.6 

l.l 

1 
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2.9 

2.3 

2.5 

^ 50p  level  means  that 

50^  of  the 

rain  events 

are  smaller 

than 

indicated. 

70%  level  means  that 

70%  of  the 

rain  events 

are  smaller 

than 

indicated,  and 
The  size  of 

so  on. 
an  event 

is  defined 

as  the  time 

for  which  rainfall 

rate  was  larger 

than  a given  value. 

multiplied  by  the  wind  velocity 

Sclontlftc  Reports  (Series  MW)  of  the  Stormy  Weather  Group 


>W-1 : Effect  of  particle  shape  and  secondary 

scattering  on  microwave  reflections  from  clouds 
and  precipitationi  by  Milton  Kerker  and  Walter 
Kitshhfeld,  March  1931  • 

MW-2:  Measurement  of  snow  parameters  by  micro- 

waves,  bv  J.S.  Marshall  and  K.L«S.  Gunn,  May 
1951. 

MW-3:  The  modification  of  rain  with  distance 

fallen,  by  E.  Caroline  Rigby  and  J.S.  Marshall, 
January  1952. 

MW-/.:  Interpretation  of  the  fluctuating  echo 

from  randomly  distributed  scatterers:  Part  I, 

by  J.S.  Marshall  and  Walter  Hitschfeld,  October 
1951. 

MW-5:  Scattering  and  absorption  of  microwaves 

by  a melting  ice  sphere,  by  M.P.  Langleben  and 
K.L.S.  Gunn,  March  1952. 

MW-6:  Interpretation  of  the  fluctuating  echo 

from  randomly  distributed  scatterers:  Part  II, 

by  P.R.  Wallace,  December  1951. 

MW-7:  The  microwave  properties  of  precipitation 

particles,  by  J.S.  Marshall,  T.W.R.  East  and 
K.L.S.  Gunn,  July  1952. 

MW-g:  Precipitation  trajectories  and  patterns, 

by  J.S.  Marshall,  M.P.  Langleben  and  E.  Caroline 
Rigby,  August  1952. 

MW-95  A theory  of  snow  crystal  habit  and  growth, 
by  J.S.  Marshall  and  M.P.  Langleben,  July  1953. 

MW-10:  The  modification  of  rain  in  showers  with 

time,  by  E.  Caroline  Rigby,  and  J.S.  Marshall, 
March  1953. 

MW-11:  A mathematical  treatment  of  random  coa- 

lescence, by  Z.A.  Melzak  and  Walter  Hitschfeld, 
March  1953. 

MW-12:  Errors  inherent  in  radar  measurement  of 

rainfall  at  attenuating  wavelengths,  by  Walter 
Hitschfeld  and  Jack  Bordan,  June  1953  • 

MW-13:  Radar  evidence  of  a generating  level 

for  snow,  by  K.L.S.  Gunn,  M.P.  Langleben,  A.S. 
Dennis  and  B.A.  Power,  July  1953. 

MW-1/.:  Initiation  of  showers  in  cumuli  by  snow, 

by  A.S.  Dennis,  July  1953. 

MW-15:  Turbulence  in  clouds  a*;  a factor  in  pre- 

cipitation, by  T.W.R.  East  and  J.S.  Marshall, 

July  1953. 

MW-16j  The  terminal  velocity  of  snow  aggregates, 
by  M.P.  Langleben,  January  195.-. 

MW-I7:  Development  during  fall  of  raindrop  size 

distributions,  by  E.  Caroline  Rigby,  K.L.S.  Gunn 
and  Walter  Hitschfeld,  January  l'^'^. 


MW-18:  The  effect  of  wind  shear  on  falling  pre- 

cipitation, by  K.L.S.  Gunn  and  J.S.  Marshall, 
December  1954. 

MW-I9:  The  convection  associated  with  release 

of  latent  heat  of  sublimation,  by  R.H.  Douglas 
and  J.S.  Marshall,  December  1954. 

MW-20:  A:  Size  distribution  generated  by  a 

random  process,  by  Walter  Hitschfeld.  B:  The 

distribution  with  size  of  aggregate  snowflakes, 
by  K.L.S.  Gunn  and  J.S.  Marshall,  September  1956. 

MW-21 : Pattern  in  the  vertical  of  snow  gene- 

ration, by  R.H.  Douglas,  K.L.S.  Gunn  and  J.S. 
Marshall,  July  1956. 

MW-22:  Precipitation  mechanisms  in  convec- 

tive clouds,  by  T.W.R.  East,  January  1956. 

MW-23i  Measurement  and  calculation  of  fluc- 
tuations in  radar  echoes  from  snow,  by 
Walter  Hitschfeld  and  A.S.  Dennis,  July  1956. 

MW-24:  The  plan  pattern  of  snow  echoes  at 

the  generating  level,  by  M.P.  Langleben, 

February  1956. 

MW-25:  A possible  role  of  hall  in  formation 
of  tornadoes,  by  Walter  Hitschfeld  and  J.S. 
Marshall,  March  1957* 

MW-26:  Growth  of  precipitation  elements  by 

sublimation  and  accretion,  by  R.H.  Douglas, 

May  1957. 

MW-27:  Studies  of  Alberta  hall  storms  1957, 
by  R.H.  Douglas  and  Walter  Hitschfeld,  May 
1958. 

MW-28:  Electronic  constant  altitude  plan 

position  Indicator  for  a weather  radar,  by 
T.W.R.  East,  .*4ovember  1958. 

MW-29:  The  motion  and  erosion  of  convective 

storms  in  severe  vertical  wind  shear,  by 
Walter  Hitschfeld,  July  19>9. 

MW-30:  Alberta  hail,  1958,  and  related 
studies.  Parts  I and  II  by  R.H.  Douglas, 

Part  III  by  R.H.D.  Barklie  and  N.R.  Gokhale, 

July  1959. 

MW-3I:  The  quantitative  display  of  radar 

weather  patterns  on  a scale  of  grey,  by 
T.H.  Legg,  June  I960. 

^n/•32:  Weather-radar  attenuation  estimates 

from  raingauge  statisltcs,  by  P.M.  Hamilton 
and  J.S.  Marshall,  January  1961. 

MW-33 : Improvements  in  weather-radar  grey 

scale,  by  F.T.  Barath,  July  1961. 
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Interim  account  of  hall  studies  - 
November  19^-0,  by  R.H.  Douglas,  J.S.  Marshall 
and  R.H.D.  Barklle,  reprinted  In  April  19o2« 

MW-35 ! Alberta  Hall  Studies,  1^61,  by  A#E. 
Carte,  RtH.  Douglas,  C.  East,  K.L.S.  Gunn, 

Walter  Hltschfeld,  J.S.  Marshall,  E.J. 

Stansbury,  December  1961. 

MW-3':;  Alberta  Hail  Studies,  1962-63,  by 
A.E.  Carte,  R.H.  Douglas,  R.C.  Srlvastava  and 
C.N.  Williams,  August  1963* 

MW-3?:  Precipitation  profiles  for  the  total 

radar  coverage,  by  P.M.  Hamilton,  September 
196/.. 

MW-3?:  Two  studies  of  convection,  by  R.C. 

Srlvastava  and  C.D.  Henry,  October  1964. 

MW-39:  Interpretation  of  the  fluctuating  echo 

from  randomly  distributed  scatters:  Part  3, 

by  Paul  L.  Smith  Jr.,  December  1964. 

MW-..0:  Facsimile  and  areal  integration  for 

weather  radar,  vols.  I and  II,  by  Marcel! 

We In,  April  I 965* 

MW-l1  : Time-dependent  characteristics  of  the 

heteorogeneous  nucleatlon  of  Ice,  by  Gabor 
Vali  and  E.J.  Stansbury,  April  1965* 

MW-/,;::  Alberta  Hail  Studies,  1<364,  by  J. 

Derome,  R.H.  Douglas,  Walter  Hltschfeld,  M. 
Stauder,  July  1965  • 

MW-43:  Attenuation  of  a parallel  beam  of  light, 

particularly  by  snow,  by  Olav  Lillesaeter,  April 
1965. 

MW-44:  Measurements  of  new  fallen  snow,  by 

K.L.S.  Gunn,  August  1965* 

M'W-L5 : Measurements  of  falling  snow,  by  K.L.S. 

Gunn  and  M.  Wein. 

MW-L6:  Experiments  on  the  nucleatlon  of  ice, 

l'^'l-63,  by  G.  Vali  and  E.J.  Stansbury,  August 
IC'15. 

MW-;.?:  Studies  of  the  formation  of  precipita- 

tion in  convection,  by  R.C.  Srlvastava  and  M. 
English,  August  1966. 

MW-L?:  Part  I of  Air  Transport  Association 

Report  "Parameters  for  Airborne  Weather  Radar", 
by  J.S.  Marshall,  C.D.  Holtz,  Marianne  Weiss, 
December  1965* 

Alberta  Hall  Studies,  1966,  by  A.J. 
Chisholm,  Marianne  English,  Walter  Hltschfeld, 
Jerry  Pell,  N'.H.  Thyer,  May  196? • 

MW-50:  Convective  overturning  and  surmer  storms, 

by  J.S.  Marshall,  January  19' 9* 

MW-51:  Measurements  of  snowfall  by  optical 

attenuation,  by  Charles  Warner,  and  K.L.S.  Gunn, 
November  F-'?. 

MW-3^;  A preliminary  microwave  attenuation 
climatology  for  the  Montreal  area  based  on  wea- 
ther radsr  data,  by  R.R.  Rogers  and  K.M.  Rao, 
Jar.uarv 


MW-53:  Showers  observed  by  stereo  cameras 

and  radar,  by  R.W.  Shaw  and  J.S.  Marshall, 
October  1972. 

MW-54:  Measurement  of  snowfall  by  radar,  by 

Paul  Carlson,  March  1968. 

MW-55:  Life  cycle  of  a suraner  storm  from  radar 

records,  by  Clifford  Holtz,  May  1968. 

MW- 56:  HARPI,  1967  - The  development  and  use 

of  a helght-azimuth-range  position  indicator,  by 

I.  1.  Zawadzki  and  E.  Ballantyne,  February  1968. 

MW-57.:  Albarta  Hall  Studies,  1967.  by  Jerry 
Pell,  P.W.  Summers  and  A.H.  Paul,  Gabor  Vali. 
December  1967* 

MW>58:  Ice  nucleatlon  relevant  to  formation 

of  hail,  by  Gabor  Vali,  December  1968. 

MW-59:  The  hall  storm  of  29  June  1967.  by  A.J. 
Chisholm,  Marianne  English  and  C.  Warner,  March 
1969. 

MW-60:  Further  studies  of  the  effects  of  pre- 

cipitation on  convection,  by  F.  Ian  Harris, 
January  1969* 

MW-61:  Quantitative  hailstorm  studies  using 

broad  vertical  beam  radar,  by  Jerry  Pell,  July 
1969. 

MW-62:  Rotating-lens  stereo  cloud  photography, 

by  R.W.  Shaw,  July  1969. 

MW-63!  Interpretation  of  the  fluctuating  echo 
from  randomly  distributed  scatterers,  Part  4. 
by  R.R.  Rogers,  September  1969. 

MW-64:  Numerical  simulation  of  large  convec- 

tive clouds,  by  Takao  Takeda,  December  1969 
(in  two  volumes). 

MW-65:  Showers  and  continuous  precipitation: 

25  July  1969.  by  Peter  P.  Zwack,  September  1971 

MW-66:  Rainfall  extreme  value  statistics 

applied  to  microwave  attenuation  climatology, 
by  J.H.S.  Bradley,  February  1970. 

MW-67:  Wind  measurements  near  Alberta  hail- 
storms, 1966-67.  by  N.H.  Thyer,  June  1970. 

MW-68:  Hall  research  at  McGill,  1956-1971. 
by  Stormy  Weather  Group,  W.F.  Hltschfeld,  ed.. 
May  1971. 

MW-69:  The  low  level  mesoscale  wind  field  of 

Alberta  hailstorms,  by  G.  Ragette,  September 

1971. 

MW-70:  Pulse  compression  for  meteorological 

radars,  by  G.i..  Austin  and  R.W.  Fetter, 

January  1972. 

MW-71 : Hail  detection  with  a polarization 

diversity  radar,  by  Brian  I..  Barge,  February 

1972. 

MW-72:  Photo  and  radar  studies  of  Alberta 

storms,  by  C.  Warner,  with  contributions  by 

J. H.  Renlck,  M.W.  Balshaw  and  R.H.  Douglas, 

April  1072. 


MW-73:  A thrt*  dimensional  numerical  model 

of  atmospheric  convection,  by  J^T.  Steiner, 

July  1972. 

MW-74:  An  Interpolation  and  extrapolation 

technique  for  precipitation  patterns,  by 
Yoondae,  0.  Ahn,  May  1972. 

MW-75:  Ice  nuclei  end  convective  storms, 

by  George  A.  Isaac,  November  1972. 

MW-76:  Measurement  and  application  of  rainfall 
autocorrelation  functions,  by  1. 1.  Zawadzkl, 
January  1973 • 

MW>77:  Rain  attenuation  studies,  by  G.  Drufuca, 

March  1973  (in  two  volumes). 


Technical  Notes 


■lWT-1;  Photography  at  the  AN/CPS-9  weather 
radar  at  Montreal  airport,  by  N.P.  Langleben 
ind  Walter  Hltschfcld,  January  1955 • 

1WT>2:  The  elevation  controller  for  CAPPI 

>peratlon  of  the  AN/CPS«9  weather  radar,  by 
r.W.R.  East.  Submitted  under  Contract  No. 
iP19(604)-1579,  October  1956. 

lWT-3:  An  optical  system  for  automatic  syn> 

.hesis  of  constant-altitude  radar  maps,  by 
I.P.  Langleben  and  W.  Denis  Gahertv,  January 
957. 

IWT-4:  On  the  measurement  of  cloud  tcmper- 

tures  from  the  ground  by  Infrared  radiation, 
y Walter  Hltschfeld,  October  I960. 

WT-5s  Two-dimensional  spectra  of  preclp- 
tatlon  patterns  by  coherent  optical  analysis, 
y 1. 1.  Zawadzkl  and  R.R.  Rogers,  September 
069. 


MWT-6s  ADA  - An  Instrument  for  real-time  dis- 
play of  microwave  attenuation  due  to  rain  by 
1. 1.  Zawadzkl  and  R.R.  Rogers,  September  1969* 

hfjT-7:  Cloud  photogrammetry,  by  J.H.  Renick 

and  R.H.  Douglas,  January  1970. 
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